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We investigate confined Levy fiights under premises of the principle of detailed balance. The mas- 
ter equation admits a transformation to Levy - Schrodinger semigroup dynamics (akin to a mapping 
of the Fokker-Planck equation into the generalized diffusion equation). We solve a stochastic target- 
ing problem for arbitrary stability index < < 2 of Levy drivers: given an invariant probability 
density function (pdf), specify the jump - type dynamics for which this pdf is a long-time asymptotic 
target. Our ("/x-targeting") method is exemplified by Cauchy family and Gaussian target pdfs. We 
solve the reverse engineering problem for so-called Levy oscillators: given a quadratic semigroup 
potential, find an asymptotic pdf for the associated master equation for arbitrary fi. 
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I. INTRODUCTION 

Many complex physical systems (like-wise non- 
physical, e.g. economic can be satisfactorily de- 
scribed in terms of the dynamics of a certain fictitious 
particle under the action of random forces (noise), orig- 
inating from its environment. Whenever we can iden- 
tify a Gaussian noise as an emergent property of the 
environment-particle coupling, the interrelated notions of 
(thermal) equilibrium, Boltzmann asymptotic probabil- 
ity density functions (pdfs) and detailed balance generi- 
cally follow. That is the case in the standard Brownian 
motion picture, based upon kinetic theory derivations, in 
the presence of (conservative) external forces. 

However, in many stochastic systems the experimental 
data show that the description based on the introduction 
of the Gaussian noise is insufhcient, since the involved 
fluctuations turn out to generate have heavy-tailed dis- 
tributions of Levy - stable type. Those distributions are 
widespread, in a broad range of systems of varied lev- 
els of complexity: physical, chemical, biological 0, 3, 
geophysical, economic That is why a deeper under- 
standing of properties of general complex systems with 
non-Gaussian noises is extremely desirable. 

For example, contrary to the case of systems with 
Gaussian fluctuations, in the context of Levy flights the 
notion of " equilibrium" , although natural under confin- 
ing conditions, has no obvious thermal connotation, see 
however Q. It is clear that any conceivable "thermal 
equilibrium" concept for non- Gaussian jump- type pro- 
cesses needs to be addressed with care and should ac- 
count for a number of precautions. In particular, an issue 
of physically motivated thermalization mechanisms for 
(confined) Levy flights has received only a residual atten- 
tion in the literature, and The main obstacle 
here may be that the source of Levy noise is interpreted 
as extrinsic to the physical system under consideration, 
with no reliable kinetic theory background, i.e. with no 
identifiable microscopic channels of an energy exchange 
with the environment. 

Levy flights are pure jump (jump-type) processes. 
Therefore, it seems useful to recall that various model 



realizations of standard jump processes (jump size is 
bounded from below and above) can be thermalized, by 
means of a locally defined scenario of an energy exchange 
with the thermostat, [ol-fTTI see also [l2|. It amounts to 
a suitable re-definition of transition rates for the jump 
process which enforces the principle of detailed balance 
to be respected by a random motion. We shall elaborate 
upon extension of this idea to Levy - stable processes, 
with a focus on the existence of asymptotic (large time 
limit) invariant pdfs, of the manifest Boltzmann form. 

Our approach is close to that used to analyze the Levy 
motion in systems with topological complexity like poly- 
mers (see, e.g. but remains distinctively different 
from a standard theory of confined Levy fiights which 
is based on the Langevin modeling. There is no known 
(additive or multiplicative) Langevin representation for 
Levy processes respecting the canonical form of detailed 
balance. 

The considered class of confined Levy fiights is well 
suited for the description of jump-type processes that are 
equilibrated (eventually, to a thermal equilibrium state) 
by a mild spatial disorder of the physical environment 
in which jumps take place. The inhomogeneity of the 
environment is quantified by turning over from the mas- 
ter equation to the affiliated semigroup dynamics. It 
is a suitable functional form of the semigroup potential 
(which we consider to be a continuous function) that al- 
lows for a unique asymptotic invariant state. That en- 
sures the existence of an asymptotic invariant pdf for the 
master equation in question. 

The structure of the paper is as follows. First we dis- 
cuss an issue of detailed balance for standard jump pro- 
cesses and next define its immediate generalization to 
Levy fiights (/i-family of Levy-stable laws with Q < ^ < 
2), Sections II and III. A mapping of the resultant mas- 
ter equation to a fractional version of the generalized 
diffusion equation follows in Section IV. For clarity of 
presentation, we make a Brownian detour in Section V 
to indicate how the semigroup framework is related to 
the standard Fokker-Planck dynamics of diffusion-type 
processes. In Section VI we describe the Levy ^ - tar- 
geting under an assumption that target pdfs are selected 
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from so-called Cauchy a - family of pdfs. For a com- 
putationally advantageous example of a = 2 and arbi- 
trary p £ (0, 2) we provide analytic formulas for the as- 
sociated semigroup potentials (they define the semigroup 
dynamics which makes the considered pdfs to be genuine 
asymptotic targets of the jump- type process). In Section 
VII the Levy targeting is considered for Gaussian tar- 
get pdfs. Section VIII presents a complete solution of 
the reverse engineering problem for the /x-family of Levy 
oscillators, corresponding to quadratic semigroup poten- 
tial. The obtained analytic formulas for asymptotic pdfs 
are depicted in Figs. 4 and 5. Not to overburden the 
paper with formal arguments, a general solution of the 
reverse engineering problem for arbitrary semigroup po- 
tential has been moved to another publication. 



Here an obvious analogue of the familiar Helmholtz free 
energy T{pt) = Y.x£K^i^)Pt{^) ~ ^(Pt) ^as been in- 
troduced, with S{pt) = —J2xeK Pt(-'')^^Pti^) being the 
Shannon entropy of the probability distribution pt{x). 
We have F{pt) > F{peq) = - InE.GX exp[-[/(.T)]. The 
relative entropy is monotonous in time and converges to 
zero, which is accompanied by a decrease of the free en- 
ergy J-{pt) to its minimal value J-{peq)- 

It is useful to mention an interesting inverse stationary 
problem of Refs. [13: El- Namely, for an arbitrary posi- 
tive probability distribution Peq{x) > on there exits 
a function U{x) such that Peq{x) is invariant under the 
jump dynamics with the transition rate kij(x,y) of the 
form (U]) . In the original formulation of Ref . , the ref- 
erence transition rate fco(a;, y) needs not to be symmetric. 



II. JUMP PROCESSES AND DETAILED 
BALANCE 

Let K he a, finite state space, with x,y € K. We con- 
sider Markovian stochastic dynamics for a finite random 
system, with transition rates k{x\y) = — > x). Given 
an initial probability distribution po{x), its time evolu- 
tion for times t > is governed by the master equation: 



dt 



Pt{x) = ^ [k{x\y)pt{y) 



k{y\x)pt{x)] 



(1) 



Given a stationary solution Pf,q (x) of the master equation, 
Peq(x) = 0. If we have 

k(.x\y)Peq{y) = k{y\x)peq{x) (2) 

one says that the condition of detailed balance is fulfilled. 

Let Peq{x) (X exp[—U{x)], where U is a suitable func- 
tion on K. (The inverse temperature /3 can be safely 
absorbed in the definition of U. As well, for clarity of 
discussion, we can set /? = 1). Accordingly: 

k{x\y)^kiy\x)eMU{y)-Uix)]. (3) 

We note that ko{x\y) = ko{y\x), in a finite state space, 
yields a uniform distribution Peq{x) = const for all x G 
K. Let us consider a simple multiplicative modification 
of a symmetric transition intensity fco(a;|y): 

'U{y)-U{x 



ko{x\y) ku{x\y) = fco(a;|y)exp 



(4) 



By inspection (simply replace k{x\y) by ku(x\y) in 
Eqs. ((H)- ©) one verifies the validity of the detailed bal- 
ance condition, with Peq{x) (x exp[— [/(a;)] as the corre- 
sponding stationary distribution. 

We assume that an equilibrium density Peq{x) > is 
unique and presume the detailed balance condition ([2|). 
([3]) to be respected. Then, the relative entropy (negative 
of the KuUback-Leibler entropy) becomes 

S{,Pt\Peq) = = J'ipt) - ^{Peq) > 0. 

(5) 



III. DETAILED BALANCE FOR LEVY 
FLIGHTS 

The above reasoning gives an immediate justification 
to the strategy adopted before in the context of Levy - 
stable processes, albeit with no explicit reference to the 
detailed balance principle, in a number of papers [Tsl - fisj . 
We also note Refs. 0, H, [l!], where "stochastic target- 
ing" and related "inverse engineering" (terms, originally 
coined in Ref. [l3l) have been exploited to this end. 

To proceed further, we recall that a characteristic func- 
tion of a random variable X completely determines a 
probability distribution of that variable. If this distri- 
bution admits a pdf p{x), we can write < exp(ipX) >= 
Jj^ p{x) exp(ipx)dx. A classification of infinitely divisi- 
ble probability laws is provided by the Lcvy-Khintchine 
formula for the exponent —F{p) of < exp(ipX) >= 
exp[-i^(p)]. 

We restrict subsequent considerations to a subclass 
of stable probability distributions with F{p) = 
with < /i < 2. The induced jump- type dynamics 
< exp(ipX4) >= eyi'p[—tF(j))] is conventionally inter- 
preted in terms of Levy fiights and quantified by means 
of a pseudo-differential (fractional) analog of the heat 
equation for corresponding pdf 

dtp = -\^\'''^ p ^ ([w^{x\y)p{y)-w^{y\x)p{x)]dy, (6) 



which has been rewritten as a master equation for a ran- 
dom system on real axis, with a pure jump dynamics. 
The jump rate w^{x\y) oc — y\^^^ is a symmetric 
function, Wfj,{x\y) = w^{y\x) akin to kQ{x\y) of the previ- 
ous subsection. We recall that the action of a fractional 
operator |A|^/^ on a function from its domain is defined 
by means of the Cauchy principal value of an involved 
integral: 

_(iA|./2 ^ np+i)M^Pi2) J dz . 

(7) 
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Mimicking the previous step (|4]), we open a possibil- 
ity of a locally controlled energy exchange with an en- 
vironment, by modifying the jump rate Wfj_{x\y) of the 
free (neither external forces nor potentials) fractional dy- 
namics to the non-symmetric form w'j{{x\y) ^ 



<(x|y) = wMv) exp([f/(y)-(7(x)]/2). With 
replacing w^(2;|?/), the master equation (j6|) ultimately 
takes a slightly discouraging form, known from a number 
of previous publications: 



dtp - ~\^\Tp - \ K {x\y)p{y)~^wl{y\x)p{x)\dy = -[exp(-[//2)] |Ar/2[cxp(C//2)p]+pexp([//2)|Ar/2 cxp(-C//2) 

^ 



The above transport equation cannot be transformed to 
any known form of the fractional Fokker-Planck dynam- 
ics, based on the standard (Levy-stable) Langevin mod- 
eling, (c.f. [l^-[2^ for literature sample). These two 
dynamical patterns of behavior are inequivalent, (3. [l6l|. 

For a suitable (to secure normalization) choice of C/(x), 
Peq(x) (X exp[— C/(x)] is a stationary solution of Eq. ([8|). 
The detailed balance principle of the form ^ , ([3]) holds 
true. 

For the record, let us mention that the free fractional 
Fokker-Plack equation ([6|) has no stationary solutions. 
Thus, the jump-type dynamics with properly modified 
jump rates clearly may give rise to confined Levy fiights. 
Their asymptotic pdfs in principle may have an arbitrary, 
not necessarily finite and/or small, number of moments. 
The reference stable laws generically have no moments of 
order higher than one. 

IV. LEVY SEMIGROUP MODELING 

The master equation ([8]) cannot be derived within the 
standard Langevin modeling of confined Levy fiights, 
16, 18|]. The latter motion scenario (with an ample 
coverage in the literature, [l9l - [2l| ) is incompatible with 
that based on the detailed balance principle ([2]), ([3]) and 
the resultant Eq. ([8]), c.f. [3,[il|. 

The form of Eq. is not handy. However, there exists 
an equivalent description of the pertinent dynamics in 
terms of a Levy-stable semigroup or a fractional (Levy-) 
Schrodinger-type equation, jl, [ij, [IE [H] . The difference 
with pure time-dependent Schrodinger equation is the 
absence of imaginary unit i before time derivative (e.g. 
inEq. ©). 

To this end let us consider the Levy-Schrodinger 
Hamiltonian operator with an external potential 

|A|"/2+V(x). (9) 

Suitable properties of V need to be assumed, so that 
— ifp is a legitimate generator of a dynamical semigroup 
exp(— iiJp) and dt^ ~ H^'^ holds true for real functions 
^{x,0) *(a;,t). 

Let us a priori select an invariant probability density 
Peq{x) = p*{x) oc exp[—U{x)] of Eq. ([S]). To make it an 



asymptotic pdf of a well defined jump-type process we 
address an issue of the existence of a suitable semigroup 
dynamics. 

Looking for stationary solutions of the affiliated semi- 
group equation dt^ = Hfi,^, we realize that if a square 
root of a positive invariant pdf (x) is asymptotically to 

1/2 

come out via the semigroup dynamics VE* — >■ p* , then the 

1/2 

resulting fractional Sturm-Liouville equation i?^p* = 
imposes a compatibility condition upon the functional 
form of V(a;), that needs to be respected. Namely, the po- 

1/2 

tcntial function and invariant pdf p, should be related 
as 

v = - ' ' ■ (10) 
p* 

The resulting semigroup dynamics provides a solution 
for the Levy stable targeting problem, with a predefined 
invariant pdf. 

Inversely, if we predefine a concrete potential function 
V{x), then the functional form of an asymptotic invari- 

1 /2 

ant pdf p* (x) (actually p* (x)) comes out from the above 
compatibility condition. We call the problem of deriva- 
tion of p* from a predefined semigrou p p otential V{x) as 
reverse engineering problem , sec Ref. [17| where this idea 
had been put forward. 

For V = V{x) bounded from below, the integral kernel 
k{y,s,x,t) = {exp[— (t — s)H]}{y,x), s < i, of the dy- 
namical semigroup exp^—tH) is positive. The semigroup 
dynamics reads: \l/(x, t) = J ^(y, s) k{y, s, x, t) dy so that 
for all < s < t we can reproduce the dynamical pattern 
of behavior, actually set by Eq. ((81), but now in terms of 
Markovian pdfs p(x, s, y, t) : 

1/2 / 

p{x,t)^pj {x)'ii{x,t)= p{y,s,x,t)p{y,s)dy, (11) 
where 

p{y, s, X, t) = k{y, s, x, t) * ^ . 

P* [y) 

An asymptotic behavior of '^{x,t) — >■ pj''(x) implies 
p{x,t) p*(x). 
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A remark is in place here. The spectral theory of frac- 
tional operators of the form Q h as received a broad 
coverage in the mathematical j23 - l28j and mathemati- 
cal physics literature [1^ [s^. An explicit functional 
form of asymptotic invariant pdfs of confined Levy flights 

1 /2 

p* (p* in the semigroup notations) is seldom accessi- 
ble, with a notable exception of those for Cauchy flights 
[2]| . Therefore it is wise to rely on accumulated data 
that are available, about the near-equilibrium behavior 
and the decay of pdfs as |a;| — >■ oo, under very general 
circumstances. Various rigorous estimates pertaining to 
the decay at infinities of the eigenfunctions, quantify the 
number of moments of the associated pdfs for different 
classes of potential functions V(a;). As well, fractional 
versions of Feynman-Kac formula determining an inte- 
gral kernel of the semigroup operator, and thence the 
transition probability which generates (by virtue of Eq. 
([TT]) ) the pdf p{x,t) dynamics consistent with Eq. 
have an ample coverage therein. 



V. BROWNIAN DETOUR 

The aim of this section is to describe the relation be- 
tween above Levy - Schrodinger semigroup framework 
and standard Fokker-Planck dynamics of diffusion-type 
processes. To make this description clear, here we put 
explicit relations, translating things from the language 
of partial differential equations (like Fokker-Planck one) 
and dealing explicitly with pdfs into the operator lan- 
guage, inherent in (both normal and fractional) quantum 
mechanics and ultimately in Levy - Schrodinger semi- 
group. 

In the theory of standard Brownian motion, the 
Langevin equation or the like (stochastic differential 
equation with the Wiener noise input) allows to infer a 
corresponding Fokker-Planck equation. This in turn can 
be transformed into a Hcrmitian (strictly-speaking, self- 
adjoint) spectral problem, (2^. Contrary to the Levy- 
stable case, for diffusion-type processes both these de- 
scriptions (e.g. semigroup and Langevin-based Fokker- 
Planck approaches) are similar descriptions of the dy- 
namics of p{x, t). 

Given the spectral solution for the operator H ~ 
— A-|- V, the integral kernel of exp{—tH) reads k{y, x, t) = 
exp(— Ej-f) $j(?/)4>*(a;). Here, the sum may be re- 
placed by an integral in case of a continuous spectrum 
and (generalized) eigenfunctions may be complex-valued. 

If we set V(a;) = identically, a purely continuous 
spectral problem arises. Then, one arrives at the familiar 
heat kernel 

k{y,x,t) = [exp{tA)]{y,x) = 

(27r)~^/^ J cxp{—p^t) exp{ip{y — x)) dp = 

(47ri)"^/2 exp 



{y - xf 



At 



which is a well-known transition probability density of 
the Wiener process (actually, upon setting t — )- (i — s)). 

When confining potentials are present, either entire 
spectrum or its part turns out to be discrete, the cor- 
responding eigenfunctions being real-valued. A standard 
example is the harmonic oscillator i.e. the Ornstein - 
Uhlenbeck process in its original stochastic version. Con- 
sider 

= (l/2)(-A + a;2-l). 

The integral kernel of cxp{—tH) is given by the classic 
Mehler formula 

k{y, X, t) = k{x, y, t) = exp{-tH){y, x) = 



7r\/r 



-2t 



cxp 



(xe * - y f 



-2t 



The normalization condition 



k{y,x,t)exp[{y^ -x')/2]dy=l 



directly employs (and defines upon setting t — > (t— s)) the 
transition probability density of the Ornstein-Uhlenbeck 
process, 

p{y,x,t) = k{y,x,t) * 

p* [y) 

with ps,{x) = Ti^^l"^ cxp(— a;2) being its (Gaussian) invari- 
ant pdf. 



VI. CAUCHY FAMILY OF PDFS AND LEVY p, 
TARGETING 



Here we describe in some detail the Levy stable 
(with stability index p) targeting strategy with the pre- 
determined one-parameter family of Cauchy target pdfs: 



p*(x) = pa{x) 



T{a) 



1 



V^r(a- 1/2) (1 



:, a> 1/2. 



(12) 

We consider functions ()12[) as asymptotic invariant pdfs 
for the stochastic jump- type process of Eq. (jS]). We wish 
to demonstrate that any /x-stable driver can be employed 
to this end. 

Instead of addressing directly Eq. (jS]), we use the semi- 
group dynamics exp(— tiJ^) generated by the fractional 
operator ([9]), i.e. the integro-differential equation 



(13) 



1 /2 

where ^'(x,i) = p{x,t)/pj {x) and V^(a;) = 

-m^^'^pi'^)/pV\Q<p<2. 

We note that the Cauchy family has been chosen 
for computational convenience only. In principle, there 



5 



is no restriction on the choice of any other target pdf 
p^{x). The quaUtative outcome will be the same as that 
provided in terms of family (fT2|) . Hereafter we call such 
general procedure " /x-targeting" . 

Let us add, as a side comment, that the Cauchy fam- 
ily of pdfs has played an important role in the previously 
mentioned search for "thermodynamic equilibria", that 
may possibly be associated with confined Levy flights, 
0-[3- It is known 0, [sj, that an exponent a in prin- 
ciple can be directly related to the thermal equilibrium 
label a oc l/fc^T. An analogous observation has been 
reported in Refs. 0, after transforming pdfs ([T^ into 
an "exponential form", which resembles Boltzmann one 
p^. (X exp(— ?7). with U{x) = Q!ln(l + x^). 

To pass over to the semigroup description we need 
to infer Vp,{x)^ given p*. This can be done analytically 
by means of the Fourier transform, specifically because 
Fourier images of functions for arbitrary a > 0.5 
exist in a closed analytical form of MacDonald functions 

The Fourier image g{k) ~ g{x)e^^^dx of a 

function g{x), when adopted to g{x) = |A|^/-^/(a;) reads 
\k\^f{k)- Fourier images of the square roots of pdfs 
read 



or (l+a\ 



7rr(a- l/2)r(a/2) 



(14) 

An explicit expression for the a-family of " - potentials" 
V^,a(x) = -{\A\''/^pi^^)/p]/^ readily follows 



2'^ 

Vp.a(a;) ^(1 + 0; 



X2F1 ( — r-^,— 



(15) 



where 2^1(0, b; c, x) is a hypergeometric function [3^. 

The expression (|15p gives the general form of the semi- 
group potentials Vp,^a{x) for arbitrary a and p. To have 
a better feeling about the properties of the function ([T^ , 
we should explore this expression for some specific values 
of parameter a. Further discussion is limited to the case 
of a = 2, i.e. 



1/2/ \ 
P2 W 



2 1 



TT 1 + X^ 



pr{k)^e-\'\. 



(16) 



1 /2 

We note that the Fourier image P2 (k) directly comes 
from the general expression (|14p . if we use -^'1/2(3;) — 
(7r/2.T)^/^ e""^. Then for all < < 2 we have 



Vfiai^) = -(1 + -^^) 2''r(l + ^)cos[(l -|-^)arctani 



(17) 



For p = I from (|17|) we recover our elder result, origi- 
nally obtained in the context of Cauchy flights, [l6l |: 



Vl,2(x) 



1 



(18) 



0.5 
go.o 

-0.5- 





/ / / 






/// 









FIG. 1; Dependence Vp..2(x) for pi terminating pdf. Figures 
near curves correspond to p values. The potentials for p = \ 
and 2 are given by Eqs. (|18p and (|19[) respectively. 



The plots of the /x-dependence of ((T7|) are reported on 
Fig.l. 

The stability index p is constrained to stay within an 
interval < < 2. The boundary value p = 2 takes us 
beyond the jump- type "territory" to continuous (Wiener 
noise) stochastic processes. It is interesting to observe 
that on the level of " //-potentials" , the transition from 
/i < 2 to /i = 2 is actually smooth. 

Analytically, recalling the fractional derivative tran- 
scription (— A)^/^ = —d^/d\x\^^ and then setting 
"blindly" p = 2 in (jlOp . we arrive at the semigroup po- 
tential for the operator H = — A + V2,2- 



V2,2(x) = Vppix) 



dx- 



1/2/ \ 
:P2 K^) 



2(3a;2 - 1) 



1/2/ \ 

P2 y.^) 



(1+X2)2 



(19) 



The notation Vfp{x) refers to the fact that this poten- 
tial appears in the semigroup (self-adjoint) version, (c.f. 
Ref. [22[) of the standard Fokker-Planck equation for 
a diffusion-type process. The same result ([T9| can be 
obtained from Eq. ([T7|) at /i = 2. 

The expression p7p permits us to expand the potential 
"1^/^,2(2^) near p = 2 to obtain 
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V,. 



>2,2 



{x) 



2(3x' - 1) 

(1+X2)2 (1+.t2)2 



[2x{x^ - 3) arctanx + {3x^ - 1) (27 - 3 + ln(l + .x^))] 



(20) 



where 7 sa 0.577216 is Euler constant. This (along with 
numerical curves from Fig.l) demonstrates the continu- 
ous transition from < 2 to /u = 2 in V^_2(2;). 



VII. GAUSSIAN TARGETING FOR LEVY 
FLIGHTS 

In the previous pubheations [1, H, [l^ we have investi- 
gated various patterns of jump-type and diffusive behav- 
ior that would produce a priori selected, basically heavy- 
tailed pdfs in the large time asymptotics. While an as- 
sociation of jump type-processes with pdfs possessing a 
finite number of moments is rather natural, an observa- 
tion of Ref. P that diffusion-type processes may as well 
admit such asymptotic pdfs, may be classified as "unnat- 
ural" . 

Here we proceed in the very same "unnatural" vein, 
asking for a Levy-stable jump-type dynamics, whose 
asymptotic pdf would have a definite Gaussian form. Let 
us select the Gaussian target pdf 



1 



p* 



aV2TT 



(21) 



whose square root pl^^{x) = f{x) = {2t:<7^) ^Z"* 
exp(—x^/Aa^) has Fourier image {p*y^^{k) = f{k) = 
(2cr2/7r)i/4 exp(-/c2cr2). That gives 



V^,G{x) = ^e^T 



1 



iFi 



1 



X 

where iFi{a, b, x) is a hypergeometric function [32[. This 
/i-family of semigroup potentials sets solution to the Levy 
stable targeting problem, if the desired target has the 
Gaussian form. 

Minor comments are necessary for a qualitative asses- 
ment of the above analytic result. The potential Vficix) 
depends on two parameters: order of fractional 
derivative p and variance a. It can be seen from Eqs. (j2ip 
and p2|) that the variance a simply alters the width of 
the potential curve and does not infiuence its shape. The 
same is true for the factor ct~^ in front of Eq. ([22]) . That 
is why in Fig. [5] we report the shape of the potential ([^^ 
in normalized variables z = x/{2<7) and ~ (T^V^G(a;). 
These universal curves are the same for any a and depend 
on the single parameter p. Note, that in these variables 
the p = 2 parabola assumes the form y2 = — 1/2. 

It is also seen from Fig. [2] that at small p the po- 
tential is around —1 (we recollect that at /i = 
Vfj.G{x) = — 1), while at larger x it has very steep growth 
like exp(z2). These steep tails flatten as p grows and 




FIG. 2: The potential (|22p in normalized variables. Figures 
near curves correspond to /i values. 



around /i = 1.5 the exponential growth of the potential 
is replaced by power-law z'^ so that at /i = 2 we have the 
correct asymptotics z'^. 



VIII. REVERSE ENGINEERING: 
ASYMPTOTIC fi - TARGETS FOR LEVY 
OSCILLATORS 

Now we pass to a detailed discussion of a particu- 
lar class of solvable examples of the reverse engineering 
problem which well illustrates the following general strat- 
egy (its full description is moved to another publication) : 
given a priori a concrete semigroup with Levy driver, in- 
fer an asymptotic pdf for the associated master equation 
®. 

Our main idea is to adopt an approach we have devel- 
oped before, [1] (see also [H,!!!]) to the Levy osciUator 
with V{x) = a;2/2 and arbitrary stability index p. 

We begin with the equation for a terminal pdf , in- 
ferred from the /x-Levy semigroup with a predefined har- 
monic potential 

v^ix)p'j^ = ^pi'^ = -|Ar/2py2^ < /i < 2. (23) 

We take Fourier images of both sides of Eq. ((23|) to obtain 
"2 1 a2/(fc) 



1 A 



2 aA:2 



(24) 
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FIG. 3: Raw solutions of Eq. (|30p (main panel) and potential 
sign k \k\^ (inset). Curves are /i-Iabeled. Arrows show the 
correspondence between potential and raw solution for given 
jx. Thick black line on the inset shows the potential for ^ = 
0.01, which has almost rectangular shape. Solution for p = 1 
corresponds to Airy function p2|) . 



its derivative at fc = 0. This is because the Eq. (|26p is of 
the second order. After that we should find the position 




FIG. 4: Normalized solutions for Fourier images of square 
roots of terminal pdfs in fc-space. Curves are ^-labeled. 



The right-hand side of Eq. 
so that 



has the form 



2|fcr/(fc). 



(25) 



The idea to solve the Eq. (|25|) for arbitrary < /i < 2 
is borrowed from Ref. jsj], where the solution for /i = 1 
had been obtained in terms of Airy functions. The 
method of Ref. [l^l is based on the consideration of ID 
Schrodinger problem with a potential being even func- 
tion of the coordinate, which implies that the correspond- 
ing eigenfunctions should be either even or odd (see e.g. 
35l l36j ) . In particular, the ground state wave function 
should be even as it does not have nodes [35|. It can 
be shown that solution f{k), defining the Fourier im- 
age of desired terminal pdf, corresponds to the ground 
state wave function of the above Schrodinger problem. 
Generalizing the method of Ref. |3J| for arbitrary fj,, we 
can show that to obtain this function for even potential 
like |fc|^ we should consider instead of ((25l) the equation 

dk-' 



2 sign/s |/s|^/(fc) or 



dk'-' 



~dW 



= 2fc^/(fc), fc>0 
-2(-fc)^/(fc), k < 0. 



(26) 



Now the scenario of obtaining the desired f{k) is as 
follows. After finding the exponentially decaying solu- 
tion of Eq. for fc > and oscillatory one at /c < 0, 
we should require the continuity of the function f{k) and 



km of the first maximum of oscillating part and shift the 
solution to the right by km so that the first maximum 
of oscillatory part is at /c = 0. Then, "chopping" the 
rest of oscillating part and refiecting the obtained piece 
about the vertical axis to obtain the even "bell-shaped" 
function. The resultant solution in the k space should be 
Fourier-inverted and squared to yield the desired termi- 
nal pdf in the x -space. 

To fulfill this scenario, we observe the following form 
of solutions of Eq. ^ for fc > and < [H. Namely, 
for fc > 




where q = {fi + 2)/2. Here J^^x) and N^^x) are Bessel 
functions and lyix) and K^ix) are modified Bessel func- 
tions, see Ref. |32| . At a; — > oo I,j{x) is exponentially 
growing function |32| while K^, [x) is exponentially decay- 
ing |32l |. On the other hand, as a; — ^ — oo the functions 
.JjAx) and Ny{x) have "needed" oscillatory asymptotics 
|33 |. This means that to have a localized pdf, we should 
leave the term with K±_ in (|27p only. Then f(k) assumes 

2g 

the following form 
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fik) 



C21J 1 (^\k\A + C22N. f^ifciOl , < 0. 



(29) 



Now we join (glue) the obtained solutions at fc = to 
secure a continuity of a function and its first derivative. 
The gluing procedure yields 



fc > 



f [cot ^J^{u) - N^{u)\ , fc < 0, 



(30) 



f{k)^C^\k\ 
where C = C12, 

--^-4^,- = ^ifcr-^ifcr*. (31) 

We note here that for the Cauchy driver, i.e. yU = 1 we 
obtain from (l30l) the result 



/(fc) = CVkKi 



'2V2, 



C 



~2¥ 



Ai 



i (2^ fc) , (32) 



known from our earlier publication 

The " raw" solutions ([50)) are plotted on the main panel 
of Fig. |3]for different values of /x. It is seen from the in- 
set that for /X — >■ (thick black line corresponding to 
/i=0.01) the potential has the shape of almost rectangu- 
lar barrier, corresponding to decaying solution (localized 
particle inside the barrier) at fc > and oscillating one 
(free particle) at fc < [3^ . We note here that for 
potentials depicted on the inset to Fig. [3] the above kind 
of solution exist only if its eigenenergy lies between the 
limiting values of a barrier at |a;| — >■ 00 [s^ H^. In this 
case the zeroth eigenenergy, which is the case for Eqs. 

and (1211) perfectly suits the problem under consid- 
eration not only for — )■ 0, where the barrier is almost 
rectangular, but also at higher /x. This explaines the fact 
that as the shape of barrier deviates from rectangular one 
at /X increase, the oscillations at fc < start to decay, the 
strongest one being at /x = 2. Also, with the growth of /x, 
the period of the oscillations lowers, the minimum being 
achieved at /i = 2 also. 

Now we find the position km of the first maximum of 
oscillating part. Equating to zero the first derivative of 
an oscillating part of (j30p we arrive at 



■KV 

cot -^J^- 



l(u) =0, 



(33) 



where v and u are defined by ([51]) . The roots of Eq. ([55)) 
can easily be obtained numerically for different /x. 

The normalization of the obtained function can be 
achieved through the condition f^{k)dk = 1 or 



2C^ 



fi{k)dk 



fi{k)dk 



1, 



(34) 



where /i and J2 denote oscillatory and decaying parts 
of Eq. ([30| respectively. Normalized solutions in the fc- 
space for different /x's are reported in Fig. 21 It is seen 
that for small fc and on the tails, the distribution func- 
tions for higher /x's run below those for smaller /x's, while 
in the intermediate fc range the situation is opposite. 

The final step of the procedure is to invert the fc-space 
solutions to the x-space and square them to obtain the 
desired terminal pdf. For general /x this procedure can 
be accomplished only numerically. 

Fig. |S] displays both the inverted functions /(fc), cor- 
responding to square roots of the inferred terminal pdfs 
(panel (a)) and those pdfs themselves (panel (b)). The 
opposite (if compared to this in the fc-space) tendency 
is seen in the x space, where the curve corresponding to 
lowest /i lies below all other curves in the small x region 
and has slowest decay. As /t grows, the central part of 
the curve rises and tails become steeper. 

Panel (c) of Fig. |5| reports a comparison between the 
shapes of functions /(fc) and f{x). The situation here is 
the same as that for the Airy function, as discussed in 
. Namely, the function in fc-space decays quicker then 
in x-space and its value at the center is larger then that 
in a; - space. We plot here the exemplary case of /x = 0.5, 
the situation for other /i is qualitatively the same. 



IX. OUTLOOK 

The next natural step in our /t-targeting procedure is 
to obtain (numerically) the dynamics of a function p{x, t) 
for Levy oscillators with different values of /t. This can 
be done both for the semigroup process ([13p and for the 
Langevin-driven one (e.g. fractional Fokker-Planck dy- 
namics). Those patterns of temporal behavior are in- 
equivalent, although both processes may terminate at 
common pdfs with a predefined decay at infinities. The 
latter pdfs may have heavy tails, but gencrically admit an 
arbitrary (finite, eventually infinite) number of moments. 

A more general problem would be that of the existence 
of terminal pdfs, after passing from the master equation 
to the (fractional) Hamiltonian dynamics ^ with an ar- 
bitrary potential V, in one, two or three spatial dimen- 
sions. 

We note that in the case of /t = 2, when frac- 
tional Hamiltonian ([9]) reduces to ordinary quantum- 
mechanical Hamiltonian operator. In the standard quan- 
tum mechanical setting (see, e.g., Refs. [s^ [s^) the 
above question is equivalent to an issue of the existence 




X, k 



FIG. 5: Panel (a) - inverted Fourier images [p*{x)]^^^; panel (b) - desired terminal pdfs at different /i (figures). Panel (c) 
compares the behavior of the functions in k and s-spaces for — 0.5. 



of bound states in a given potential. The quantum me- 
chanical language appears because we can convert the 
parabolic equation of the Fokker-Planck type to the gen- 
eralized Schrodinger equation. 

The wave function of a bound state should be localized 
to ensure a normalization of its squared expression, i.e. 
the corresponding stationary pdf of the Fokker-Planck 
equation. It is known (see, e.g., Ref. [s^) that in ID 
case the bound state exist in the potential well U{x) of 
not only finite but an infinitesimal depth. The only re- 
striction is that the integral U{x)dx should exist. 
The latter condition is equivalent to the requirement that 
U{x) should have the same asymptotics at infinities and 
potential zero point U(±oo) = 0. In the 2D case, when 
the potential U = U{x,y), the situation is similar to that 
in ID one, while in 3D ([/ = U{x,y^z)) the situation is 
to some extent opposite - if the potential well is not suffi- 



ciently deep (see Ref. HI] for details), the particle cannot 
be "captured", so that bound state does not exist. Con- 
fining potentials in 3D, where bound states exist, form 
the so-called Kato class of potentials. 

The presence of fractional derivatives with < ^ < 2 
alters the picture both in ID (2D) and in 3D. In ID they 
definitely "spoil" the bound states. It is not only that 
the pdfs (if in existence) may have heavier tails if com- 
pared to the conventional (/j, = 2) case. The pdfs in 
question may not exist at all, if a normalizability of the 
bound state is lost. In 3D and in equations with frac- 
tional derivatives there may typically be no normalizable 
bound states (and thus terminal pdfs), except for a care- 
fully selected (Kato)-subclass of conceivable potentials. 

Some peculiarities pertaining to the (non)-existence 
of invariant pdfs in the case of Levy drivers (Langevin- 
driven fractional dynamics) were discussed for ID case 
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in Ref. [2j]. We have encountered the same problem in 
connection with the Cauchy family of pdfs [1, Q , see also 



Ref. [37[ for a discussion of so-called infinite covariant 
densities. 
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